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O ; Abstract 

Ch \ A study of general properties of the narrow-width approximation (NWA) with polarization/spin 

decorrelation is presented. We prove for sufficiently inclusive differential rates of arbitrary resonant 

^. 

CN . decay or scattering processes with an on-shell intermediate state decaying via a cubic or quartic 

(— I \ vertex that decorrelation effects vanish and the NWA is of order T. Its accuracy is then determined 
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numerically for all resonant 3-body decays involving scalars, spin-2 fermions or vector bosons. We 






specialize the general results to MSSM benchmark scenarios. Significant off-shell corrections can 

occur - similar in size to QCD corrections. We qualify the configurations in which a combined 
(N 
^ ' consideration is advisable. For this purpose, we also investigate process-independent methods to 
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improve the NWA. 
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I. INTRODUCTION 

Theoretical arguments and experimental observations indicate that new particles or inter- 
actions play an important role at the TeV scale, which will soon become directly accessible 
at the Large Hadron Collider (LHC). In the near future we can therefore anticipate ground- 
breaking discoveries that reveal physics beyond the Standard Model (BSM) jl|. Theoretically 
appealing extensions of the Standard Model often feature numerous additional interacting 
heavy particles. BSM phenomenology is hence characterized by particle production and 
cascade decays, which lead to many-particle final states and scattering amplitudes with 
complex resonance structure. In order to extract the additional Lagrangian parameters of 
an extended theory from collider data, theoretical predictions are required that match the 
experimental accuracies. 

The treatment of unstable particles in quantum field theory is complicated by the absence 
of asymptotic final states for particles with finite mean lifetime. That a consistent theory 
respecting unitarity, renormalizability and causality can be formulated based on external 
stable states only has been shown in Ref. 2j. That work, however, does not discuss how to 
carry out fixed-order perturbative calculations involving unstable particles. Subsequently, 
several finite-width schemes have been discussed in the literature that in principle facilitate 
cross section calculations for arbitrary processes at tree [3] and loop level J], where the 
full fixed-order scattering amplitude is taken into account.^ If the width F of an unstable 
particle is much smaller than its mass M, amplitude contributions that feature the corre- 
sponding resonant intermediate state are generically enhanced by a factor of order M/T. 
A theoretically consistent method to extract this typically dominant part of the amplitude 
is provided by the narrow-width approximation (NWA) 6|. The NWA allows to neglect 
nonresonant as well as nonfactorizable amplitude contributions, thus leading to significant 
simplifications for calculations in extensions of the Standard Model or calculations of higher 
order corrections. When the NWA is applied in such calculations the unstable particle is 
effectively restricted to on-shell states. We note that the NWA is hence implicitly applied 
whenever branching ratios are extracted from reaction rates. Recently, it has been observed 
that the NWA can be unreliable in relevant circumstances, namely with decays where a 



^ A method for the inclusion of finite- width effects in event generators is discussed in Ref. 



daughter mass approaches the parent mass or when parton distribution functions are con- 
volved with a resonant hard scattering process [7| ■ A more detailed study of the accuracy 
of the NWA is thus well motivated. 

In Sec. ini we describe the NWA and clarify the treatment of polarization and spin 
correlations. We then consider correlation effects for decay and scattering rates in the 
special case of an on-shell intermediate state, which allows us to prove that the NWA error 
is of order F. To the best of our knowledge no proof for this result has previously been given 
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in the literature. The polarization/spin decorrelation has been considered in Ref. 
authors find a result similar to that derived in Sec. IIIBl but obtain it by performing the 
little-group integral for the resonant particle and assuming that the subamplitudes transform 
under some suitable representation of the little group. In our more direct derivation we will 
explicitly argue for the transformation of the subamplitudes and see that the full little-group 
integral is not generally required. 

When applied in perturbative calculations, the uncertainty of the NWA is commonly 
estimated as between ^ T/M/S and ~ ST/M, i.e. of order T/M in the physics sense. 
In Sec. mil we test this assumption systematically for all resonant 3-body decays involving 
scalars, spin-i fermions or vector bosons. We then specialize the generic results to benchmark 
scenarios in the Minimal Supersymmetric Standard Model (MSSM). Our findings in Sec. lIIII 
naturally raise the question of process-independent improvements of the NWA, which we 
investigate in Sec. llVi We close with a summary in Sec. |Vl 

II. NWA PROPERTIES 

A. NWA definition 

Formally, the NWA factorization into production of the unstable particle and its subse- 
quent decay is obtained from the full cross section formula by factorizing the phase space 
and integrating out the Breit-Wigner resulting from the squared denominator of the reso- 
nant particle's propagator. With P denoting the sum of the incoming particles' momenta 
and q and pi denoting the momenta of the intermediate resonant particle and the final state 



particles, respectively, the phase-space factorization of the n-particle phase-space element 

n n ,3 -► 

into a j-particle phase-space element for the production process, an {n — j + l)-particle 
phase-space element for the decay and an additional g^-integration is given by 

d^ri{P;Pi,---,Pn) = d<^j(P;pi,...,pj_i,q) x -— x rf$„_j+i(g;pj, . . . ,p„) . (1) 

ZTT 

Denoting the resonant particle's mass and width by M and F, respectively, the squared 
propagator denominator is given by D{q^) := 1/ [(g^ — M^)^ + M^F^]. It is integrated out 
via 



^ ^ [ dq^S{q' - M^)ar{q') = ^4^, 
2MF J ^ ) rVH 1 2MF ' 

where gminmax ^^6 ^^6 kinematic bounds for q^ and (Tr{q^) denotes the residual differential 
cross section (dependencies on quantities other than q^ have been suppressed). We note 
that the second and third expression are asymptotically equal for F — > 0. For finite F, 
Transformation ([2]) constitutes the core NWA. In the first step one assumes that adding the 
contributions from q^ regions that are not kinematically accessible has a negligible impact, 
which is the case if M^ is not close to (?minmax ^s measured by F. In the second step we 
assume that the g^-dependence of the remaining integrand crr(g^), i.e. the residual squared 
amplitude and phase-space factors, is weak compared to that of the strongly peaked Breit- 
Wigner. 

While the core NWA already provides considerable simplification, it does not facilitate 
complete factorization for nonscalar particles, since the production and decay process are 
linked by correlation effects. By additionally neglecting polarization/spin correlations be- 
tween the production and decay parts we obtain the standard NWA. More specifically, 
the corresponding numerator of the squared matrix element is replaced for unstable vector 
bosons and fermions, i.e. 
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FIG. 1: Resonant decay or scattering process kinematics with total incoming momentum P. 

with vector boson propagator in unitary gauge and suppressing the argument q for the 
polarization vectors and spinors. The resulting factorization of the cross section into ap x BR 
constitutes the standard NWA factorization. Two sources of potential problems thus arise: 
the factorization of the squared amplitude neglects correlation effects, and off-shell effects 
are ignored when replacing the g^-integration with a constant factor. We address the former 
issue in the next section. 

B. Polarization/spin decorrelation 



In this section we show that there are no polarization/spin correlation effects for total 
or sufficiently inclusive differential rates^ of arbitrary resonant decay or scattering processes 
with an on-shell intermediate state decaying via a cubic or quartic vertex. We consider 
processes with an unstable vector boson and note that a similar argument applies for unstable 
fermions. The squared amplitude for a process with resonant massive vector boson, which 
we assume to be on mass-shell, is given by 
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where the denominator of the squared propagator is suppressed. By decoupling the polariza- 
tion sum and inserting a polarization average, i.e. applying Transformation ([3]), one obtains 
the decorrelated squared matrix element 
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^ For a cubic (quartic) decay vertex, 9i (6*0, ^o and ^i) has to be integrated out. The angular variables are 
defined below Eq. ([5]). 



As first step, we consider the process of Fig. [H where the resonant particle is produced in 
an arbitrary process with total incoming momentum P and decays via a 2-body decay. To 
obtain the total rate for a specific process, the matrix elements have to be integrated over 
the full phase space. Here, we consider the total rate after application of Transformation ([2]), 
which - up to overall factors - is given by 

n ,3 ^ n 
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The standard NWA is now obtained by decorrelating the squared amplitude, i.e. by replacing 



I Air- 1 with |A^rP. To determine the effect of this procedure we consider the difference 
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(4) 
Exploiting that H defined in Expression (j4]) is Lorentz invariant, we choose the Gottfried- 
Jackson frame |9i], in which q = {Eq,0)'^, to evaluate it. In this frame our assumption 
g^ = M^ implies q = (M, 0)"^. By integrating out p2 and substituting d^pi = |pip d\pi\ dQi 
we find 



H = p-^^^^^ 6{E, - E, - E,) jdn, (\Mr\'-\M, 



■ (5) 

{p2 = -pi,Eq = M} 



We choose dQi as innermost integration, for which \pi\ and p^, . . . ,pn are invariables. Since 
the production and decay amplitudes A^pC* and M.'^tu are Lorentz invariant, Aif^ and A4'^ 
are Lorentz 4-vectors. The spatial component of the decay (production) 4-vector has to be 
proportional to pi (a function of Ps, ■ ■ ■ ,Pn) since no independent other 4-vectors exist in 
the decay (production) part, i.e. with J^p = {J^^, A4p)'^ and Aid = {A4'^,Aid)'^ we have 
A4d oc pi and A^p fixed. ^ We therefore choose dQi = dcos9i dcpi with 9i = l{pi,Aip). 
Furthermore, in the Gottfried- Jackson frame the polarization sum —Qfj^u + q^qu/M"^ is given 



The proportionality factor for Add can depend on pi only through |pi|, which is fixed. 



by diag(0, I3) and for the squared matrix elements we then get 
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\Mr\^ = \Mp ■ Md\^ = \Mp\^\Md\^ cos^ 01 and \Mr\^ = -\Mp\^\Mci\^ . 
Evaluating I defined in Expression ([5]), we thus obtain 

I = \Mp\^\Md\' [ dpi [ dcosOi fcos^9i-^^=0, 

which proves that polarization decorrelation does not introduce an error for on-shell inter- 
mediate states in resonant processes of the type shown in Fig. [TJ 

We now consider the case where the unstable particle decays via a quartic vertex. More 
specifically, we consider the process of Fig. [1] with an additional particle with outgoing 
momentum po attached to the decay vertex connecting q, pi and p2. Again, we consider a 
total or sufficiently inclusive differential rate after application of Transformation ([2]) and, 
similar to Eq. ([5]), arrive at 

d^Po d^pi 
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as test for correlation effects. Now, |po|, IpiI, ^1 and M.p which, as before, is a function 
of p3,...,Pn are fixed for the inner integrations over 6q, (J)q and 0i. We parameterized 
the po-integration as d^po = \po\'^ d\po\ dcosOodcpo, where the orientation of po is nieasured 
relative to that of Aip, i.e. 6q = l{po, Aip). The pi-integration was parameterized as d^pi = 
|PiP d\pi\ (icos6'i dcpi, oriented such that 9i = l{pi,Po). ri, the unit vector in the direction of 
Pi, is obtained from ro, the unit vector in the direction of po, by first rotating by 61 around 
the axis defined by 0o = (— sin 00, cos 00, 0)'^ and then by (pi around ro: 

Pi = \Pi\ Qiro, (pi) Q{(pQ, 9i) ro , (6) 

where Q{x, a)ij = — J2k ^ijk^k sin a + {6ij — XiXj) cos a + XiXj is the matrix representing a 
rotation by the angle a around the axis defined by the unit vector x = (xi, X2,Xs)^. Since 
the moduli of po and pi as well as their relative orientation are fixed for the integrations in 
J, the decay matrix element can be decomposed as 

Md = C0P0 + ci pi + C2P0X pi , (7) 



where we have exploited that the integration is performed in the Gottfried- Jackson frame and 
Aid can thus not depend on ps, ■ ■ ■ ,Pn- Note also that p2 = —po — Pi is not an independent 
3- vector and hence does not appear in the decomposition. The coefficients Cq, Ci and C2 can 
only depend on \pq\, \pi\ or Pq -pi (or equivalently ^1). They are, however, independent of the 
variables of the J- integration 6'o, 0o a^nd 0i. Using the explicit expressions of Eqs. ^ and 
([7]), J can be evaluated, and we find J = 0. We conclude that also no error is introduced 
by decorrelation when the intermediate particle decays via a quartic vertex. 

We now generalize the proof to an on-shell intermediate state decaying via a cubic or 
quartic vertex in an arbitrary resonant process, where the decay products of the unstable 
particle may decay further. Suppose the particle with momentum pi decays into m particles 
with momenta pn, . . . , pim- We map the kinematics of the extended process to the kinematic 
configuration of Fig. [1] and factorize the entire phase-space element by applying Eq. ([1]) to 
the intermediate state pi 

d^P; ...) = d^{P;pi,p2, ...,pn)^ d^puPn, ■ ■ • ,Pim) • 

By first integrating d^{P;pi,p2, . . . ,Pn), "we construct the situation that we dealt with be- 
fore, the only difference being that in the decomposition of the matrix element Aid we now 
have additional terms involving pu, i = 1, . . . ,m. However, in the crucial inner integrations 
of d^{P;pi,p2, ■ ■ ■ ,Pn) the relative orientations of these momenta are fixed and thus they 
transform like pi under rotations. The absence of on-shell correlation effects can therefore 
be shown in the same way as above. Generalization to the case where multiple decay prod- 
ucts decay further or where correlation effects for multiple on-shell intermediate states are 
considered is straightforward. Note that this extension again holds not only for total rates, 
but also for sufficiently inclusive differential rates. 

C. Order of the NWA error 

We can now prove that the relative error 

j3 ToFS — TnwA /qn, 

^ = r ^^^ 

J- NWA 



is of 0{r) for the decay and scattering rates considered in Sec. IIIBI ^ A priori this is not 
clear. The error as function of F could, for example, behave like V^ with < k < 1. In this 
case, a series expansion around F = would be invalid. We consider the off-shell decay rate 



^D(g^F)F.(g^), 

^min 

where kinematic hmits define the integration bounds ^mmmax ^^^ ^^ have exphcitly noted 
the dependence of the squared propagator denominator on F. Fr(g^) denotes the residual, F- 
independent integrand.^ More specifically, F^ is the squared combination of the production 
matrix element, the numerator of the propagator and the decay matrix element, integrated 
over d$p and d^^ and includes overall factors like l/(2vP^). The decay rate in NWA is 
given by Fnwa = rp(M^) x Vd{M'^)/V where Fp and F^ are the production rate and partial 
decay rate, respectively. We can now write: 



ToFs-rNWA = / -— £'(g^F)F 



'- '^D{q\V) (F.(g2) - V ,{M')2MT ,{M^)) - aiV) ^p^^")^'^^^") 
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rymax fin 

/ -^Z}(g^r)(F,(g2)-F,(M2))-a(F) 



with a := 1 - 2MF /t'"" £'(g^ V)dq^/{2'n). Note that a = for g^j^ ^^^ = ±00. As shown 
in Sec. IIIBl decorrelation effects vanish and Fr(M^) factorizes, which has been used in (*). 



For the relative deviation R we get 
ToFS ~ Tnwa 
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F 



NWA 
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/ 9max rtrt^ 
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F.(g^) - F.(M2) , , 



Fp(M2)F,(M2) 



Assume now that Tr{q'^) is twice continuously differentiable in the kinematically allowed 
phase-space region and that Fp(g^)Frf(g^) = <^ g^ G {(Jmin; Q'max}- Then, the function 

Tr{q'^) 



^ We illustrate the proof using total decay rates, but note that it also applies to sufficiently inclusive 

differential decay and scattering rates as considered in Sec. Iff Bl 
^ Since interference effects with nonresonant diagrams are not resonance enhanced, we do not consider them 

here. 



is twice continuously differentiable with respect to q^ for all q^ G [gmim Q'max] ^'^'^ ^^ ^ 
(^min) Q'max)) wliicli we assuuie in the following, since we consider resonant intermediate states. 
By second order Taylor expansion we then get 

9M{q') = 9m{M^) + 9'm{M^) ■ {q' - M^) + \g'Ux) ■ {q' - M^ 

with X G [o'min' 9'max]- Notc that X depends on q^. As continuous functions on the compact 
set [^min) 9max]) Q'm ^-^d g'j^ are then bounded and we can find Li > and L2 > such that 



9M{q)\ < Li and |^M(r)l < ^2 W e [qii^,q, 



(10) 



We can now further evaluate \R\ to obtain 

\R\ < 
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By employing the Taylor expansions 



^°Sy^-^ = logy + C'(a:') 
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X 2 k 



for small positive x to expand the log and arctan functions for small F, we get 
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Lastly, we evaluate and expand a: 
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and obtain as final result: 
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where L^ and L2 are F-independent constants with mass dimension —1 and —3, respectively. 
We have thus shown that R is of 0(r) if M is in the interior of the kinematically allowed 
region. RjY is then finite in the limit F -^ 0. Note that this is a stronger statement 
than the asymptotic equality of Fqfs and Fnwa for F ^ 0, which immediately follows from 
the asymptotic equality of the off-shell and core NWA expressions (see Transformation ([2])) 
and the absence of on-shell correlation effects. R = 0(r) does, however, not guarantee 
that \R\ ~ T/M as suggested by the scales that occur in the unstable particle propagator. 
The leading terms in F on the right-hand side of Ineq. ( fTTl) . while F-independent, may 
nevertheless become arbitrarily large as M approaches the kinematic limits. This can be 
seen directly by inspecting the factors that are explicitly (iJminmax ~ ^^) dependent, but is 
also due to the fact that Li and L2 may become large as M^ approaches g^j^ or g^ax- 

III. NWA ACCURACY FOR RESONANT THREE-BODY DECAYS 

In Sec. [TTl we have already determined several necessary conditions for a small relative 
NWA error \R\: The standard NWA will only be accurate if polarization/spin correlation ef- 
fects can be neglected. Furthermore, Transformation ([2]) shows that Ko'minmax)^''^ ~ ^\ ~ ^ 
is necessary to effectively eliminate the dependence on the g^-integration bounds. In the 
next section we systematically explore NWA deviations for the process class that allows ap- 
plication of the NWA while featuring minimal complexity in order to identify other accuracy 
limiting factors. 

A. Model-independent analysis 

We have systematically probed the NWA accuracy for resonant 3-body decays involving 
scalars, spin-^ fermions or vector bosons. By inserting these particles into the topology of 
Fig. [2] one obtains 48 generic decay processes. A particular process is selected by giving type 
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FIG. 2: Resonant 3-body decay kinematics. 
(4-monientum, mass) for each particle: 

Tj{Pi,Mj) -^Tiipumi), T{q,M) and T{q,M) -^ T2{p2,m2), Ts{p^,m3), 

where type can be scalar (S), fermion (F) or vector boson (V). The width of the intermediate 
particle with momentum q is F. We have created a custom, model-independent program 
that allows to calculate total decay rates with and without NWA for all 48 processes. The 
program has been validated using SMadGraph 10| and SDECAY Uj.^ Using this program, we 
scanned the parameter space for large values of |-R| by varying masses, width and couplings. 
Note that coupling constants typically cancel in R with the exception of the relative strength 
of the chiral components of SFF and VFF vertices. This dependence has been taken into 
account in the scan. Qualitatively, we find that the NWA error will exceed order T/M for 
mass configurations in an extended vicinity of kinematic bounds, but not for sufficiently 
central configurations. A detailed description and discussion of the scan and its results can 
be found in Ref. [l3j. Here, we choose the process that features the largest deviations, SSS- 
SSV, to explain the mechanism that can lead to large deviations even when the g^-integral 
is not noticeably cut off by kinematic bounds. The largest deviations occur when all but 
one final state mass are small. We therefore set mi = m2 = 0, but keep rris ^ 0. The 
g^-integrand for this process is then given by 
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(t 



Mj 



rris 



2\2 



(r - ml) 



(12) 



q- J \ ml ) (g2 - M^)^ + M'^Y'^ ' 

where the 1st- and 2nd-stage decay phase-space elements contribute the first and second 
factor, respectively, and the 2nd-stage decay matrix element gives the third factor. When 



^ Unitary gauge has been applied tliroughout. Tlie employed Feynman rules have been double-checked 
using Ref. [12 1. 
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FIG. 3: Resonant 1^3 decay SSS-SSV (see main text): The graph displays the g^-dependence 
of the Breit-Wigner that is integrated out in the NWA (dashed) and of the complete integrand of 
Eq. (I12p (solid) for m^ = M — 3T. The contour plot shows R' := |-R|/(r/M), i.e. the magnitude 
of the relative NWA error in units of T/M, as function of m^ and M with rrii = ?ti-2 = 0. The 
width-mass ratio F/M is 0.01. 

ml ^ M^ the second and third factor effect a strong deformation of the Breit-Wigner 
shape, which, together with the resulting large NWA deviations, is displayed in Fig. [3l The 
deviation grows with increasing power of the deforming factors. When M approaches the 
lower kinematic bound, |-R| is sensitive to the type of the 2nd-stage decay, which determines 
the power of the factor that deforms the Breit-Wigner peak. While this factor enhances 
the Breit-Wigner tail, the factor of the Ist-stage decay suppresses it. And vice versa for 
the upper bound. Generally, we find stronger effects for SSV, VSV, FFV, VVV and SVV 
than for FSF, SFF, VFF, VSS and SSS vertices. As mentioned above, coupling parameters 
do not affect the relative NWA error if they factorize. However, for processes with SFF or 
VFF vertices, a strong dependence on the relative strength of the chiral components can 
exist. In Fig. HI we demonstrate this for the FFV- VFF decay. We conclude this section with 
a caveat: Fofs has been calculated including the full range of invariant masses g^ for the 
intermediate resonant state. In a model-independent analysis we could, however, not take 
into account interference with amplitude contributions with different resonance structure 
or with nonresonant diagrams. In gauge theories, such contributions can in principle be 
important. In the next section we will consider examples in specific models. 
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FIG. 4: Dependence of the magnitude of the relative NWA error R' (see Fig. [3]) on the strength of 
the chiral components of the Ist-stage decay vertex ^'^{aPi + PPr) for process FFV-VFF at the 
parameter point M/Mj = 0.68, mi/Mj = 0.3 and m2 = rn^ = with T /M = 0.01. Color code as 
in Fig. El 

B. MSSM analysis 

We now specialize our model-independent analysis to benchmark scenarios in the MSSM. 
We determine for each 1^3 decay^ in the MSSM and each benchmark parameter point 
of the Snowmass Points and Slopes (SPS) [14] whether the resulting process is resonant, in 
which case the NWA accuracy R of Eq. ([8]) is calculated. The SPS mass spectra and decay 
widths have been generated using SOFTSUSY [15|] and SDECAY, respectively. As example, 
we select the most thoroughly studied point SPS la (defined by the parameter set Mq = 
100 GeV, Ml/2 = 250 GeV, Aq = -100 GeV, tan/3 = 10, /i > 0) and display in Table [T] 
resonant 3-body decay modes with \R\ > SF/M for which the distance between M and 
kinematic bounds is larger than 5 widths. Note that for bi -^ W^ti -^ W^bxt the NWA 
error R is negative. All other processes in Ta ble [II have positive R. A complete listing of 
results with \R\ > 5T/M can be found in Ref. 13|. 

The tabulated results are intended as guideline to alert the NWA user of potential large 
errors. In such cases a more detailed on- /off-shell comparison including all amplitude con- 



^ We consider final states with stable particles and up to one unstable particle. This is of interest, since 
subsequent decays do not generally mitigate the deviations [7|. 
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TABLE I: Resonant 3-body decay modes in the MSSM at SPS la, for which R' , i.e. the mag- 
nitude of the relative NWA error in units of T/M (see Fig. [3]) is larger than 5, and A := 
min (M/ — ttt-i — M, M — m2 — m^) /T, i.e. the minimal distance between M and the kinematic 
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tributions is called for. We exemplify this by discussing two specific cases. The first process 
is ^ ~^ Xidd^ which proceeds via intermediate dL,R states (see Fig. [5]). At SPS la, one has 
Mg = 607.7 GeV, Mj^ = 545.2 GeV, Mj^ = 568.4 GeV and M^o = 96.7 GeV, which implies 
that all intermediate states can be on-shell. As seen in Tabled one finds R ~ lOF/M for 
the graph with intermediate d^.^ In contrast, the graph with intermediate dpt is not affected 
by an unexpectedly large NWA error {R ^ 2 T/M) . This can be traced to the fact that d^ 



Note the equality of the total decay rates with intermediate cIl.r and d^ ^ state due to the C invariance 
of|X|2. 
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X? ~9 ^--^-^^X'^ 

dh dn 



d d d d 



X? ~9 ^-.-^^X? 

dl d\ 



FIG. 5: Complete set of Feynman graphs for the MSSM 3-body decay g -^ Xi^d. 

has a much higher decay probabiUty than d^ due to the chiral structure of the chargino and 
neutrahno couphngs. Consequently, the total width over mass ratio is about 1% for di, but 
only about 0.05% for dn, i.e. in the latter case KfJminmax)^^^ — M\ > 200T and hence the 
absence of large NWA errors. When also taking into consideration the much larger resonant 
enhancement for the dji in comparison with the di mediated process, one can expect that 
the former process in NWA provides a good approximation to the full g — > x^'^rf rate. When 
calculated in NWA, the decay rate of the di mediated process is indeed approximately 100 
times smaller than the decay rate of the dji mediated process. Under these circumstances 
interference effects are likely to provide the leading correction to the NWA estimate. In- 
terference between amplitudes with intermediate states related by charge conjugation is 
Breit-Wigner suppressed, because almost always either the invariant mass of Xid or Xid will 
be close to M, i.e. only one of the propagators will be resonant. Furthermore, interference 
between amplitudes with di and dR intermediate states is highly suppressed, since the over- 
lap of the corresponding resonances is negligible: M^ — Mj ^ Fj + Tj . We computed 
the total decay rate for g -^ x\dd by Monte-Carlo integration of the SMadGraph-generated 
complete matrix element. Comparing this off-shell result to the NWA result for the dn 
contribution, we find a relative deviation of 1.2%, which confirms our expectation. 

As second example, we consider the process xt ~^ Xi'^d, which proceeds via intermediate 
iy+, ul and dl states (see Fig.E]). At SPS la, one has M~+ = 181.7 GeV, M^o = 96.7 GeV, 
Mw = 79.8 GeV, M^^ = 561.1 GeV and Mj^ = 568.4 GeV, which implies that only 
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FIG. 6: Complete set of Feynman graphs for the MSSM 3-body decay Xi ~^ Xi^d. 

the intermediate W~^ state can be on-shell. For the Ist-stage decay Xi ~^ Xi^^ with^ 
BR = 7.5% the minimal distance to the kinematic bounds is Mj — mi — M = 2. OF, which is 
not large enough to rule out an unexpectedly large NWA error. In fact, Table [I] shows that 
R = 13% = 5.2 F/M for xt ~^ Xi^^ ~^ Xi'^d. Primarily via interference the nonresonant 
Ul and d*j^ contributions slightly reduce the NWA error to i? = 11%, where the total decay 
rate for xt ^ Xi'^d with complete matrix elements is compared to the xt ~^ Xi^~^ ~^ Xi'^d 
rate in NWA. Since the Ist-decay stage is not affected by QCD corrections, this NWA error 
is particularly significant. 

IV. PROCESS-INDEPENDENT NWA IMPROVEMENT 

In this section we propose a modified NWA prescription that preserves the simplifying 
NWA features and improves the behavior close to kinematic bounds by taking into account 
phase-space properties. The standard NWA is obtained by replacing the distribution ip : 
f *~~^ I '2^^iQ^)fi^^)i "which - acting on the residual differential cross section ar - gives the 
total cross section^°, by a rescaled Dirac distribution, i.e. 

^(^r) = j ^Diq')ariq') ^ ^NWA(a.) = ^j dq^5{q^ - M')ariq') , (13) 

which is motivated by the observation that the Breit-Wigner shape suppresses the contribu- 
tions with g^ 7^ M^. That is, for small F, ip essentially eliminates all contributions of f{q^) 
except for the on-shell part, and in the limit F — > 0, ip^ar) and i^NWAicTr) are asymptotically 
equal. 

However, for finite F the g^-dependence of the phase-space factors and residual matrix 
elements can cause a significant deformation of the shape of D{q^)ar{q^) relative to the 



At SPSla' [l3] the branching ratio is 1.3%. 

Apphcation to cross sections and decay rates is analogous. 
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Breit-Wigner shape D{q^), resulting in a considerable shift of the maximum position and 
maximal value. The effect is particularly strong when M is close to the kinematic bounds, 
where threshold-type factors suppress the resonant contribution and shift the maximum. 
Focusing on the impact of the g^-dependence of the process-independent phase-space factors, 
we write ilj{(Jr) = J ^D{q'^)FS{q'^)ar{q^), where PS(g^) denotes the integrand factor arising 
from the phase-space element. Even though the g^-dependence of PS(g^) deforms the Breit- 
Wigner, the shape of D(g^)PS(g^) is also strongly peaked and thus suggests a NWA-inspired 
approximation, which takes into account the shift of the maximum position caused by PS(g^). 
More specifically, we propose to substitute the mass of the resonance M with an effective 
mass Meflf in ip^wAi'^r)- In analogy to M^ being the maximum position of the Breit-Wigner, 
Mgg is given by the position of the maximum of D(g^)PS(g^): 



^NWA(crr) -^ ^PSiNWA(o-r) = ^,. ^ / dq^6{q'^ - M^g)crr{q^ 

ZiWeffi J 

1/2 

with Meff := I argmaxD(g^)PS(g^) 



(14) 



The thus defined effective mass only exploits kinematic information and is hence universal 
for the class of processes with identical phase-space properties. ^^ Due to limr^o^eff = M, 
the deviation of the effective mass from the physical mass is negligible unless M is close 
to the kinematic bounds. Therefore, the proposed phase-space improved narrow-width ap- 
proximation (PSINWA) does not result in significant deviations in cases where the standard 
NWA gives (9(r/M)-accurate results. On the other hand, when M approaches a kinematic 
bound the behavior is improved: Since the distance between Mes and the bound stays finite 
the PSINWA result does not vanish. The PSINWA error is therefore bounded in contrast 
to the diverging standard NWA error. 

To exemplify this method, we consider the scalar resonant 3-body decay SSS-SSS (see 



Sec. nil Ap . With (3{m^M) := a/1 — m? /M'^ and kinematic conventions as in Fig. [2], the 
3-particle phase-space element is given by 



dq^ 



^^ Since M^g can be obtained by numerical, one-dimensional maximization of D{q^yPS{q^) with q^ ~ M^ 
as suitable initial value, the computational complexity of its determination is negligible. In theories with 
mass relations it may be necessary to adjust other parameters to maintain theoretical properties like gauge 
invariance. 

18 



R 


\ 
\ 


8 




6 


\ 
\ 




\ 


4 


\ 




^^^ \ 


2 


^ \Co^. 






-2 


M/(m2 +m3) 




o.i^ 



M/{Mi -mi) 



FIG. 7: Comparison of the phase-space improved (sohd) and standard (dashed) NWA error for 
a scalar resonant 3-body decay, i.e. process SSS-SSS of Sec. IIIIAl when the resonance mass M 
approaches the kinematic bounds m2 + rris (left) and M/ — mi (right) (see Fig. [2]). Displayed is 
the relative approximation error in units of resonance width/mass, i.e. R = (FoFs/rpsilNWA ~ 
l)/(r/M), for T/M = 0.05, mi = ms = and m^/Mi = 0.2. 
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M^Q is thus obtained via maximization of 






PS(g2)D(g2) oc 



In Fig. U\ we compare the error of the phase-space improved and standard NWA when the 
kinematic bounds are approached. It shows that the effective-mass method mitigates the 
error |i?| to less than SF/M for M arbitrarily close to the kinematic bounds. The divergence 
of the standard NWA error in this region is also clearly visible. With increasing distance to 
the kinematic bounds, on the other hand, the PSINWA and standard NWA results converge 
and arebothof 0(F/M). 

The PSINWA is general in the sense that the effective mass does not depend on the process 
matrix elements, and as effective-mass method it preserves the on-shell simplifications of 
the standard NWA. The former strength is, however, also a potential shortcoming, since 
for some processes Breit-Wigner deformation is also caused by the g^-dependence of the 
residual matrix elements, which may have to be taken into account to achieve the desired 
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approximation accuracy. As long as the deformed Breit-Wigner shape is still strongly peaked, 
the proposed method can be extended: a process-specific effective mass can be calculated 
by maximizing /(g^) := -D(g^)PS(g^)|A^(g^)P (and possibly also treating the width as free 
parameter ).^^ For cases where the Breit-Wigner deformation is so severe that a modified 
mass is not sufficient to obtain the desired accuracy, it may be possible to successfully apply 
the process-specific method proposed in Ref. [l7|. 

V. SUMMARY 

We studied the general properties of the NWA with polarization/spin decorrelation. After 
defining the NWA and clarifying the treatment of polarization and spin correlations, we 
proved for sufficiently inclusive rates of arbitrary resonant decay and scattering processes 
with an on-shell intermediate state decaying via a cubic or quartic vertex that decorrelation 
effects vanish and the NWA is of order F. When applied in perturbative calculations, 
the NWA uncertainty is commonly estimated as between ^ T/M/S and ~ 3F/M. We 
tested this assumption by systematically determining the NWA accuracy numerically for 
all resonant 3-body decays involving scalars, spin-| fermions or vector bosons. We found 
that the approximation error will exceed the (9(F/M) estimate for mass configurations in 
an extended vicinity of segment kinematic bounds. This is due to a significant distortion of 
the Breit-Wigner peak and tail, which is effected by the g^-dependence of the phase-space 
elements and the residual matrix elements. While factorizing coupling parameters do not 
affect the relative NWA error, for processes with SFF or VFF vertices a strong dependence 
on the relative strength of the chiral couplings can exist. We specialized the general results 
to MSSM benchmark scenarios and presented results for SPS la. We found that significant 
off-shell corrections can occur - similar in size to QCD corrections. To simplify a combined 
treatment we proposed a modified, process-independent approximation that exhibits an 
improved accuracy close to kinematic boundaries while preserving the simplifying power of 
the standard NWA. 
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For more complicated phase spaces and/or matrix elements the function to maximize, namely f{q'^), may 

ff 



not be available in analytical form. The computation of Af?™ will then be more expensive 
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